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—The  gross  spectral  character  and  duration  of  the  acoustic  wave 
signature  of  a  nuclear  explosion  near  the  ocean  surface  is  estimated 
by  constructing  theoretical  pressure-time  histories,  using  models 
for  the  explosion  and  wave  propagation.  The  explosion  is  assumed  to 
have  a  yield  of  1  KT  and  the  nominal  range  is  6600  km.  The  fre¬ 
quencies  of  interest  are  rather  low,  50  Hz  and  less,  so  absorption 
is  small  and  the  estimates  of  spectral  character  and  duration 
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20.  ABSTRACT  (continued) 

essentially  decouple.  The  spectrum  depends  almost  entirely  on 
the  source  characteristics  and  the  duration  is  controlled  by 
characteristics  of  the  travel  path. 

The  airblast-induced  pressure  loading  on  the  ocean  surface 
dominates  the  source,  with  the  acoustic  waves  from  direct  couplinc 
into  the  water  being  relatively  small.  At  large  distances,  the 
spectrum  for  1  KT  peaks  near  20  Hz  and  is  band-limited  between 
5  and  50  Hz.  For  different  energy  yields  these  frequencies  scale 
with  the  cube-root  of  the  yield.  Different  assumptions  about  the 
(laterally  homogeneous)  oceanic  sound  profile  lead  to  differing 
estimates  for  the  signal  duration.  Values  of  20  to  60  seconds 
seem  most  reasonable.. 

A- 


_ UNCLASSIFIED _ 

SICUMTY  CL  Altl  FI  CATION  OF  THIS  PA0Cr*>>«<  Dm*  Cf»*rM) 


AFTAC  Project  Authorization  No.  VT/0712/B/PMP 
ARPA  Order  No.  2551,  Program  Code  No.  6H189 
Effective  Date  of  Contract  November  17, 1978 
Contract  Expiration  Date:  November  15, 1981 
Amount  of  Contract  $1315.437 
Contract  No.  F08606-79-G0008 
Principal  Investigator  and  Phone  No. 

Dr.  J.  Theodore  Cherry,  (714)  4530060 
Project  Scientist  and  Phone  No. 

Mr.  Brian  W.  Barker,  (202)  325-7581 


This  research  was  supported  by  the  Advanced  Research  Projects 
Agency  of  the  Department  of  Defense  and  was  monitored  by 
AFTAC/VSC,  Patrick  Air  Force  Base,  Florida  32925,  under  Contract  No. 
F086Q679-C0008. 

The  views  and  conclusions  contained  in  this  document  are  those  of 
the  authors  and  should  not  be  interpreted  as  necessarily  representing 
the  official  policies,  either  expressed  or  implied,  of  the  Advanced 
Research  Projects  Agency,  the  Air  Force  Technical  Applications 
Center,  or  the  U3.  Government 


Bv _ _  _ _ 

DistriT;',4.  o  ' 

_ Avail..:  ill4-  Codes 

jAvaii  ,  r.ii/or 
Dist  Special 


TABLE  OF  CONTENTS 


Section 


I.  INTRODUCTION  .  1 

1.1  INTRODUCTION .  1 

II.  SYNTHETIC  ACOUSTIC  RECORDS  .  4 

2.1  INTRODUCTION .  4 

2.2  AIRBLAST  SOURCE  .  4 

2.3  DIRECT  COUPLING  .  6 

2.4  CALCULATION  OF  ACOUSTIC  WAVES  FROM  AN 

EXPLOSION  .  6 

2.5  THE  GREEN'S  FUNCTIONS  -  ACOUSTIC 

WAVES  WITH  WKBJ  RAY  THEORY . 10 

2.6  THE  SPECTRAL  SIGNATURE  AT  THE  RECEIVER.  15 

2.7  ESTIMATES  OF  SIGNAL  DURATION . 17 

2.8  CONCLUSIONS . 23 

REFERENCES . 26 


ill 


SYSTEMS.  SCIENCE  AND  SOFTWAm 


LIST  OP  ILLUSTRATIONS 

Figure  Page 

1.  Schematic  description  of  the  construction 

of  synthetic  acoustic  wave  records .  2 

2.  Pressure  loading  for  the  airblast  from  a 

1  KT  near-surface  explosion  .  5 

3.  Spatial  distribution  of  the  pressure  applied 

at  the  surface  to  represent  the  direct 
coupling  of  a  1  KT  explosion  into  the  water  .  7 

4.  The  explosion  source  spectra .  9 

5.  Rays  included  in  the  WKBJ  ray  solution.  ...  13 

6.  The  spectrum  of  the  pressure  time  history 


of  a  single  ray  arrival  at  a  SOFAR  axis 
receiver  station . 16 

7.  Sound  velocity  profiles  .  18 

8.  The  Green's  functions  for  the  velocity 

profiles  of  Figure  7 . 19 

9.  The  dependence  of  the  Green's  functions 

on  range  for  two  models . 21 


10.  Green's  function  variations  for  different 

bottom  depths  and  near  surface  velocities  .  .  22 


iv 


SYSTEMS.  SCIENCE  AND  SOFTWARE 


I.  INTRODUCTION 


1.1  INTRODUCTION 

Nuclear  explosions  detonated  on  or  near  the  surface  of 
the  ocean  are  expected  to  excite  acoustic  waves  that  may  be 
detected  at  large  distances.  In  this  report  we  attempt  to 
estimate  the  gross  spectral  character  and  duration  of  the 
waves  associated  with  such  an  event,  assumed  to  have  a  yield 
of  1  KT. 

To  estimate  the  character  of  the  acoustic  signature  of 
the  explosion,  we  construct  theoretical  pressure-time  his¬ 
tories,  using  models  for  the  explosion  and  associated  acoustic 
wave  propagation.  The  assumptions  of  the  model  are  outlined 
in  some  detail  in  Section  II  and  are  shown  in  schematic  form 
in  Figure  1.  Basically,  we  convolve  an  elastic  representation 
for  the  source  with  a  Green's  function  representing  the  propa¬ 
gation.  The  source  may  include  two  parts.  The  first  is  the 
airblast  loading  applied  to  the  surface  and  the  second  is  the 
direct  coupling  of  thermal  energy  into  the  water.  When  the 
explosion  is  assumed  to  be  above  a  certain  height,  estimated 
to  be  about  2.3  meters  for  1  KT,  the  direct  coupling  contri¬ 
bution  is  small.  Thus,  our  calculations  include  the  cases 
when  the  explosion  is  on  or  above  the  surface,  though  we  have 
not  attempted  to  account  for  burst  height  effects  that  may  be¬ 
come  important  when  the  explosion  is  more  than  10  meters  above 
the  surface.  The  spectra  of  both  the  airblast  and  direct 
coupled  portions  of  the  source  are  quite  similar  in  character. 
They  each  are  nearly  flat  at  low  frequency  and  have  a  corner 
frequency  that  occurs  between  20  and  30  hertz. 

For  the  propagation  to  ranges  of  3000  kilometers  and 
larger,  we  use  a  single  sound  velocity  profile  and  a  wave  propa¬ 
gation  technique  adapted  from  seismology  called  WKBJ  seismogram 
method  (Chapmen,  1978,  Brown,  et  al. ,  1980).  The  assumption  of  a 
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SOURCE : 

AIRBLAST 


Plus 


//n''v 


DIRECT  COUPLING  FOR  1  KT 
NEAR  SURFACE  EVENT 


PROPAGATION: 


WKBJ  RAY  THEORY  IN  A 
LATERALLY  HOMOGENEOUS 
OCEAN 


Figure  1.  Schematic  description  of  the  construction  of 
synthetic  acoustic  wave  records.  The  source, 
represented  by  two  pressure  distributions 
applied  to  the  ocean  surface,  is  convolved 
with  a  Green ’ s  function  representing  the 
propagation. 
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laterally  homogeneous  ocean  with  a  planar  bottom  is,  of  course, 
a  rather  idealized  model  for  the  actual  situation. 

At  the  frequencies  of  interest  (<  50  Hz) ,  absorption 
in  the  ocean  is  very  small.  Therefore,  the  spectral  character 
of  the  source  should  be  preserved  at  great  distances.  On  the 
other  hand,  the  duration  of  the  acoustic  signal  is  almost 
entirely  dependent  on  characteristics  of  the  travel  path. 
Reasonable  estimates  for  the  path  parameters  suggest  that  the 
duration  is  between  20  and  60  seconds  for  a  range  of  6600  km. 

We  are  not  aware  of  empirical  data  describing  the 
character  of  the  acoustic  signatures  of  nuclear  explosions 
detonated  over  the  deep  ocean.  However,  we  have  seen  sono¬ 
grams  made  from  hydrophone  recordings  of  two  French  explosions 
detonated  high  above  a  Pacific  atoll.  The  explosion  yields 
were  of  the  order  of  half  a  megaton  and  the  burst  heights  were 
about  500  meters.  The  recordings  were  from  a  range  of  about 
6000  km.  The  sonogram  displays  of  the  hydrophone  recordings 
of  these  two  events  showed  signals  sharply  limited  in  frequency 
and  time.  The  time  duration  was  about  13  seconds  for  one  of 
the  events  and  8  seconds  for  the  other,  which  had  smaller  yield. 
The  spectra  seemed  to  have  a  distinct  corner  frequency  at  20 
to  22  hertz. 

We  are  not  sure  of  the  analog  between  the  high  yield, 
large  height-of-burst  atoll  explosions  and  low  yield  near¬ 
surface  explosions  in  the  deep  ocean,  so  relevance  of  the  re¬ 
cordings  of  the  French  tests  is  unclear.  Nevertheless,  the 
characteristics  we  are  attempting  to  bound  for  the  prototype 
event  are  the  signal  duration  and  frequency  content.  We  find 
that  the  acoustic  signal  at  large  distances  from  a  1  KT  event 
should  have  a  spectrum  that  peaks  near  20  hertz.  Nearly  all 
the  signal  energy  is  between  5  and  50  hertz.  The  duration  is 
estimated  to  be  20  to  60  seconds.  These  estimates  do  not  seem 
too  inconsistent  with  the  data  from  the  French  tests. 
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II.  SYNTHETIC  ACOUSTIC  RECORDS 


2.1  INTRODUCTION 

In  this  section  we  describe  the  assumptions  of  our 
analysis  and  summarize  the  results.  Our  conclusions  are 
summarized  in  Section  2.8. 


2.2  AIRBLAST  SOURCE 

The  most  important  component  of  a  near  surface  source 
is  the  airblast  loading  of  the  surface.  This  was  represented 
by  the  Air  Force  Weapons  Laboratory  (AFWL)  nuclear  blast 
standard  developed  by  Needham,  et  al.  (1975) .  This  model 
provides  the  pressure  loading  at  the  surface  for  a  one- 
kiloton,  free-air  detonation  at  sea  level.  The  pressure 
loading  function  is  as  indicated  in  Figure  2.  Within  a 
radius  of  5  meters  the  peak  pressure  is  70  kbar.  The  peak 
then  decays  proportional  to  rjl  when  rp  >  5  meters.  After 
the  shock  front  passes  by,  the  pressure  decays  rapidly  to 

Pin* 

In  applying  the  standard  airblast  load,  large  deforma¬ 
tion  of  the  water  is  ignored  and  acoustic  wave  theory  is 
assumed  to  be  valid.  At  each  point  on  the  surface  a  pressure, 
P(r,t)  calculated  from  subroutines  given  by  Needham,  et  al. 
(1975),  is  applied  and  the  resulting  acoustic  waves  are  com¬ 
puted,  as  wili  be  subsequently  described. 

The  standard  airblast  should  be  a  good  representation 
of  the  source  for  burst  heights  between  2.3  and  10  meters. 

The  radius  of  the  "burnout"  sphere  is  about  2.3  meters  for 
1  KT,  so  below  this  height  direct  coupling  of  thermal  energy 
into  the  water  should  be  included.  Above  10  meters  there  may 
be  height-of-burst  effects  that  should  be  accounted  for  in 
the  P (r, t) . 
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We  do  not  know  what  might  be  the  effect  of  ignoring  the 
large  "crater"  that  will  be  formed  in  the  water  by  the  near 
surface  explosion.  A  far  more  detailed  and  expensive  calcula¬ 
tion  would  be  required  to  estimate  it. 


2.3  DIRECT  COUPLING 

For  near  surface  explosions,  substantial  amounts  of 
energy  may  be  directly  coupled  into  the  ocean.  Details  of 
this  coupling  will  depend  on  the  specific  device  design,  but 
we  can  make  some  good  estimates  of  the  gross  features.  Numer¬ 
ous  previous  studies  for  the  Defense  Nuclear  Agency  (DNA)  (e.g, 
Allen,  et  al . ,  1973;  Allen  and  Baker,  1974)  indicate  that  five 
to  ten  percent  of  the  device  energy  is  coupled  to  soil  or 
water  for  a  surface  burst.  Therefore,  we  completed  a  non¬ 
linear,  hydrodynamic  one-dimensional  (plane  strain)  calculation 
of  a  100  ton  nuclear  explosion  in  water.  The  calculation  was 
continued  into  the  linear  response  region,  which  was  estimated 
to  begin  at  a  radius  of  15  meters.  A  pressure  pulse  of  the 
form  shown  in  Figure  3  was  then  chosen  as  an  "equivalent  elastic 
source."  This  pressure  pulse,  when  applied  over  10  msec,  gives 
the  same  total  impulse  as  the  detailed  nonlinear  calculation. 

The  radius  of  20  meters  is  the  radius  where  the  airblast 
loading  dominates  the  direct  coupling  contribution  predicted 
by  this  calculation.  Note  that  this  is  not  much  different 
from  the  15  meter  elastic  radius. 


2.4  CALCULATION  OF  ACOUSTIC  WAVES  FROM  AN  EXPLOSION 

The  airblast  and  direct  coupling  portions  of  the  source 
are  both  represented  by  pressure-time  histories,  P(r,t) 
applied  to  the  ocean  surface.  Each  point  on  the  surface  is 
then  a  source  of  acoustic  waves  and  the  total  solution  is 
obtained  by  convolving  with  the  appropriate  Green's  functions. 
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In  this  case,  the  Fourier  transformed  pressure  at  the  receiver, 

A 

Pr  (ui) ,  is  given  by 

Pr(w)  »  iui  nR  G(ui)  F  («j)  ,  (1) 


where 


cD  is  the  sound  velocity  at  the  source  and  p  is  the  ray  para- 
meter.  The  effective  Green's  function,  iajnRG(aj),  is  the  spec¬ 
tral  pressure  at  the  receiver  due  to  a  vertical  point  load 

A 

applied  at  the  center  of  the  source  region.  The  F  (ui)  is  the 
total  force  applied  by  the  source. 


F(oj) 


(r,oj) 


JQ(upr) 


rdr 


(2) 


where  P(r,uj)  is  the  Fourier  transform  of  the  applied  pressure, 
which  is  assumed  to  be  negligible  for  r  >  b.  All  calculations 
were  done  with  a  single  ray  parameter  (p  =  0.667  sec/km).  That 
is,  we  ignore  small  variations  in  F(lj)  for  ray  parameters  (take¬ 
off  angles)  slightly  different  from  this  one.  For  each  ray 
this  formulation  is  valid  for  distances  that  are  large  com¬ 
pared  to  the  source  dimension,  and  for  the  time  window  of  in¬ 
terest,  all  rays  have  essentially  the  same  geometrical  ray  para¬ 
meter. 

A 

In  Figure  4  we  plot  the  source  spectra,  F(id),  associated 
with  the  airblast  and  direct  coupling  contributions.  For  an 
ocean  path  G(w)  is  expected  to  be  a  weak  function  of  u>  (though 
modulated  by  interference  from  many  separate  arrivals) ,  and 
the  pressure  spectrum  at  the  receiver  should  look  roughly  like 
oj  times  F  (w) . 
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An  important  feature  of  the  F(w)  is  that  they  both  are 
low  frequency  sources  that  roll-off  sharply  above  25  or  30 
Hertz.  The  airblast  source  is  about  three  times  as  large  as 
the  direct  coupling  contribution.  Thus,  as  far  as  acoustic 
waves  are  concerned,  there  is  little  to  distinguish  between 
sources  on  and  above  the  surface. 

The  amplitude  of  the  spectra  in  Figure  4  are  controlled 
by  the  total  impulse  for  each  source  component.  In  view  of 
our  simplifying  assumptions,  the  relative  values  are  certainly 
only  rough  estimates  for  the  impulse  amplitudes.  However, 
they  are  probably  correct  in  indicating  that  the  airblast  and 
direct  coupling  contributions  are  the  same  order,  with  the 
direct  coupling  contribution  being  somewhat  smaller. 

The  corner  frequencies  in  the  F  (ui)  are  controlled  by 
the  source  dimensions.  Again,  our  calculations  give  rough 
estimates  for  this  quantity,  though  errors  of  as  much  as  a 
factor  of  two  would  be  surprising. 

2.5  THE  GREEN'S  FUNCTIONS  —  ACOUSTIC  WAVES  WITH  WKBJ 

RAY  THEORY 

For  the  acoustic  wave  propagation  in  the  ocean  sound 
profile  we  used  a  method  called  WKBJ  seismogram  method.  This 
method  was  developed  by  C.  H.  Chapman  (1978)  for  computing 
synthetic  seismograms  in  earth  models.  (A  more  accessible 
description  of  WKBJ  ray  theory  is  given  by  Aki  and  Richards, 
1980,  Section  9.4).  M.  J.  Brown  of  Scripps  Institute  of  Ocean¬ 
ography  has  adapted  this  theory  for  acoustic  wave  propagation 
in  the  ocean.  Brown,  et  al . ,  (1980)  apply  the  technique  to 
model  220  Hz  acoustic  wave  propagation  at  a  range  of  900  kilo¬ 
meters  in  the  northwest  Atlantic.  The  results  indicate  that 
this  technique  gives  some  improvement  over  results  from  geo¬ 
metric  ray  theory. 
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Given  the  pressure,  Pr (w) ,  at  the  receiver  due  to  a 
concentrated  pressure  PQ  at  the  source,  the  Green's  function 
is  computed  from 

pc{z)  2P  (id) 

G(u>)  -  - ^ — - -  ,  (3) 


p  is  the  density  at  the  receiver  and  c(z)  is  the  sound  velocity 
which  depends  on  the  depth,  z.  The  Pr(oj)  is  the  Fourier  trans¬ 
form  of 


Pr(t)» 


4II2(2  r  )1/2 


Im 


d 

at 


A(t) 


E 

j 


E 

t-0 


P1/2  R.. (p,Z) 

| <3© j  (p,r,z)/dp| 


(4) 


where  R^(p,z)  is  the  product  of  reflection  coefficients,  which 
will  be  further  described  below,  and 


Alt) 


H (t)  .  i  H(-t) 


(5) 


9j(p,r,z)  »  pr  +  t  j  (p,  r ,  z)  .  (6) 

In  these  equations,  r  is  the  range  of  the  receiver.  The  travel 
time  function  is 


tj(p,z)  = 


I 

ray path 


(7) 
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> 


» 


The  first  sum  in  Equation  (4)  is  over  all  j  rays  that  reach 
the  receiver.  The  second  sum  is  over  all  p  that  satisfy 
t  -  6 j  (p,r,z) . 

Geometrical  ray  theory  assumes  a  simple  form  for  the 
contributions  due  to  ray  parameters  near  the  geometrical  ray 
parameter  in  equation  (4) ,  and  replaces  the  second  sum  with 
a  simple  time  function.  The  WKBJ  seismogram  includes  these 
contributions  more  rigorously.  It  is  usually  considered  to 
be  an  excellent  approximation  for  periods  much  less  than  the 
total  travel  time. 

The  Rj{p,z)  is  given  by 


Rj  (p*z) 


n 


(8) 


where  is  the  product  of  the  reflection  coefficients 

along  the  ray.  These  are  -i  for  each  turn  within  the  ocean, 

-1  for  a  reflection  at  the  surface  and  £  <  1  for  each  reflec¬ 
tion  at  the  ocean  bottom. 

The  types  of  rays  included  in  the  WKBJ  theory  are 
indicated  schematically  in  Figure  5.  The  parameter  £  accounts 
for  leakage  of  energy  into  the  bottom  sediments  on  each  bottom 
reflection.  Estimates  for  this  bottom  loss  are  available  in 
the  literature,  though  at  frequencies  somewhat  higher  than 
those  of  interest  here.  Using  a  plot  in  Section  5.8  of  Urick 
(1975)  we  estimated  that  0.977  should  be  on  upper  bound  for 
5. 

Absorption  in  the  water  is  very  small  at  the  fre¬ 
quencies  (<  50  Hz)  of  interest.  It  was  ignored.  Also  ignored 
is  any  loss  from  reflection  at  the  surface. 

Shown  at  the  right  in  Figure  5  is  a  typical  sound 
velocity  profile.  Indicated  on  this  profile  is  the  point  at 
depth  which  has  the  same  velocity  as  the  velocity  at  the 
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BOTTOM  REFLECTIONS 


) 


DIFFRACTION 


f 


.gure  5.  Rays  included  in  the  WKBJ  ray  solution.  A  typical 
sound  velocity  profile  is  shown  at  the  right.  The 
dot  indicates  the  point  at  which  the  velocity  is 
the  sane  as  at  the  surface. 
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surface.  All  rays  must  turn  below  this  depth.  As  currently 
implemented,  this  theory  is  unable  to  account  for  waves  trapped 
in  the  region  where  velocities  are  lower  than  at  the  surface; 
that  is,  SOFAR  channel  waves. 

The  most  serious  limitation  of  the  theory  is  that  it 
fails  to  account  for  scattering  due  to  lateral  changes  in  the 
ocean  bottom  and  in  the  sound  velocity  profile.  Such  scatter¬ 
ing  would  introduce  energy  into  the  low  velocity  channel,  where 
it  would  be  trapped  and  travel  quite  efficiently.  The  theory 
is  being  modified  to  account  for  such  phenomena. 

In  summary.  Equation  (1)  indicates  that  the  acoustic 
wave  signature  of  the  event  is  obtained  from  the  derivative  of 
the  convolution  of  the  source  function  (Figure  4)  with  the 
Green's  function,  G(ui),  would  be  essentially  constant  over  the 
frequency  band,  though  with  some  scalloping  due  to  interference 
between  arrivals,  then 

Pr(w)  a  i  oj  F  (in)  (9) 

gives  a  smoothed  estimate  for  the  spectrum  at  the  receiver. 

The  WKBJ  theory  is  a  finite  wavelength  theory  that 
gives  higher  order  corrections  to  geometric  ray  theory.  How¬ 
ever,  numerical  tests  indicate  that  the  only  gross  differences 
between  G(ui)  computed  with  our  WKBJ  seismogram  program  and  that 
expected  from  geometric  ray  theory  were  due  to  the  low  frequency 
band  width,  which  is  determined  by  discretization  of  the  time 
series  and  the  smoothing  required  to  prevent  aliasing.  Thus, 
as  far  as  the  important  characteristics  of  the  answers  are 
concerned,  we  may  as  well  have  used  geometric  ray  theory. 

Our  final  estimates  for  the  acoustic  waves  from  the 
near  surface  explosion  were  made  in  the  following  way.  The 
spectrum  was  computed  from  Equation  (9).  That  is,  the  travel 
path  plays  no  role  in  filtering  the  spectrum,  other  than 
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through  the  differentiation  that  comes  from  the  force-pres¬ 
sure  relationship.  This  is  the  result  we  would  get  with  geo¬ 
metric  ray  theory.  The  duration  of  the  recorded  signal  is 
estimated  by  the  WKBJ  ray  theory  program.  However,  as  we 
will  describe  in  Section  2.7,  we  will  make  some  crude  adjust¬ 
ments  to  these  estimates  in  an  attempt  to  account  for  scat¬ 
tering  effects  not  present  in  the  model. 

2.6.  THE  SPECTRAL  SIGNATURE  AT  THE  RECEIVER 

As  we  explained  in  the  previous  section,  for  frequencies 
less  than  60  Ht  or  so,  the  oceanic  acoustic  waves  are  propa¬ 
gated  with  little  frequency  distortion.  Thus,  (9)  gives  an 
estimate  fcx  the  (smoothed)  spectrum  of  the  signal  recorded 
at  large  dif. ranees  from  a  nuclear  explosion  on  or  near  the 
ocean  surface.  That  is,  we  simply  multiply  the  spectra  in 
Figure  4  by  iai. 

In  Figure  6  we  plot  our  estimates  for  the  spectra  of 
the  pressure  pulse  at  a  hydrophone  station.  The  plots  with 
and  without  the  direct  coupling  contribution  are  shown.  As 
expected,  it  makes  little  difference. 

The  other  important  characteristic  of  the  spectra  in 
Figure  6  is  that  they  are  band  limited,  having  a  peak  near 
22  Hz  and  being  down  a  factor  of  3  from  this  peak  at  about 
4  and  50  Hz. 

One  factor  we  have  not  included  in  these  calculations 
is  the  hydrophone  response.  At  low  frequencies,  the  response 
is  roughly  proportional  to  u>  (LTC  T.  Jerrick,  AFTAC,  oersonal 
communication) ,  which  makes  the  spectrum  fall-off  even  more 
sharply  on  this  side.  The  hydrophone  response  on  the  high 
frequency  side  of  the  peak  is  not  known  to  us. 

We  should  also  point  out  that  the  spectra  cube-root 
scale  with  yield.  That  is,  if  the  yield  were  twice  the  1  KT 
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Figure  6.  The  spectrum  of  the  pressure  time  history 
of  a  single  ray  arrival  at  a  SOFAR  axis 
receiver  station.  The  complete  signal 
spectrum,  which  is  a  sum  of  many  single  ray 
spectra,  would  have  the  same  basic  shape. 
The  amplitude  units  are  arbitrary,  depen¬ 
ding  on  the  amplitude  of  the  Green's  func¬ 
tion  for  the  travel  path. 
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value  assumed,  the  spectral  peak  would  move  from  22  Hz  to 
17.5  Hz  and  the  region  with  amplitudes  within  33  percent  of 
the  peak  moves  from  4-50  Hertz  to  3  -  40  Hertz. 

2.7  ESTIMATES  OF  SIGNAL  DURATION 

To  estimate  the  signal  duration,  we  computed  the  Green's 
functions  for  several  typical  oceanic  profiles.  The  Green's 
function  is  the  pressure  at  the  receiver  due  to  a  concentrated 
force  on  the  surface  at  the  source.  This  was  computed  using 
the  WKBJ  seismogram  technique  described  in  the  previous 
section. 

The  ocean  profiles  were  taken  from  compilations  of 
sound  velocity  data  for  the  south  Atlantic  Ocean.  Four  pro¬ 
files  were  chosen  that  represent  extremes  during  the  early 
Fall  in  this  part  of  the  ocean.  These  profiles  are  shown  in 
Figure  7. 

The  Green's  functions  were  computed  at  ranges  from 
3600  to  6600  kilometers.  At  6600  km  a  typical  calculation 
includes  rays  with  90  to  140  turning  points.  There  are  a 
similar  number  of  bottom  bounces.  As  mentioned  before,  we 
assumed  the  reflection  coefficient  at  the  bottom  to  be  £  * 
0.977.  Then  for  100  bounces  the  total  attenuation  is  0.01. 
Therefore,  the  bottom  bounces  make  almost  no  contribution  to 
the  record  at  such  large  ranges. 

In  Figure  8  we  plot  the  Green's  functions  at  6600 
kilometers  for  the  four  profiles  in  Figure  6,  with  the  bottom 
depth  assumed  to  be  5000  meters.  The  records  are  dominated 
by  the  delta  function  contributions  from  the  turning  rays. 
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Figure  7.  Sound  velocity  profiles. 
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The  bottom  bounces  arriving  earlier  and  later  are  much  smaller. 
There  are  also  some  diffraction  effects,  but  they  are  quite 
small. 

The  duration  of  the  Green's  functions  was  estimated 
from  the  window  within  which  the  amplitudes  are  30  percent  of 
the  maximum  and  greater.  The  measurements  were  done  by  eye, 
with  no  attempt  to  be  precise.  For  example,  the  duration  for 
the  "average"  model  in  Figure  8  was  taken  to  be  9  seconds. 

Similar  calculations  were  done  for  several  parameter 
variations.  In  particular,  we  varied  the  following: 

1.  Range  -  The  duration  was  found  to  be  directly 
proportional  to  the  range,  as  should  be  ex¬ 
pected.  The  Green's  functions  are  shown  in 
Figure  9. 

2.  Bottom  Depth  -  Shallower  depths  give  a  smaller 
window  for  the  turning  rays  (Figure  5) .  There¬ 
fore,  there  are  fewer  such  rays  and  a  shorter 
duration.  Deeper  depths  give  a  larger  window 
and  commensurately  larger  duration.  Several 
Green's  functions  illustrating  this  are  given 
in  Figure  10 . 

3.  Surface  Velocity  -  The  velocity  near  the  surface 
is  probably  quite  variable  along  the  path.  This 
velocity  was  decreased  slightly  (from  1525  to 
1520  m/sec  in  the  top  10  meters)  in  the  "average" 
path  model  (A) .  The  resulting  model  is  called 
"A  MOD."  The  main  effect  is  again  to  increase 
the  depth  window  for  turning  rays  and  thus  to 
increase  the  duration.  Some  examples  are  shown 
in  Figure  10. 
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As  we  have  mentioned,  the  failure  to  account  for  energy 
that  is  scattered  into  the  SOFAR  channel  is  the  main  weakness 
of  the  theory.  The  onset  time  is  probably  reasonably  accurate, 
if  the  profile  represents  a  good  average  for  the  entire  path. 

As  a  crude  estimate  for  the  extended  duration  due  to  late 
arriving  energy  which  propagates  in  the  SOFAR  channel,  we  com¬ 
puted  R/c0,  the  arrival  time  for  a  wave  at  the  axis.  An 
"extended  duration"  was  then  taken  to  be  the  time  between  the 
onset  time,  measured  on  the  Green's  function  plots,  and  R/cQ. 

In  Table  1  we  summarize  our  duration  estimates.  If  we 
use  the  results  directly  from  the  model,  the  duration  varies 
from  4  to  38  seconds.  Extending  the  duration  by  assuming 
arrivals  at  the  SOFAR  axis  velocity,  the  duration  estimates 
are  less  variable,  ranging  from  22  to  55  seconds.  The  latter 
estimates  are  much  more  realistic. 


2.8  CONCLUSIONS 

Our  main  conclusions  about  the  acoustic  waves  at  large 
ranges  from  a  1  KT  near-surface  explosion  are  as  follows: 

•  For  the  low  frequency  (<  50  Hertz)  signals  of 
interest,  estimates  for  the  spectral  content 
and  duration  may  be  decoupled.  The  former  is 
almost  entirely  controlled  by  the  source  and 
the  latter  by  the  sound  velocity  profile  along 
the  travel  path. 

•  The  direct  coupling  of  energy  into  the  water 
seems  relatively  unimportant  compared  to  the 
airblast  loading. 

•  The  acoustic  signal  spectrum  is  peaked  at 
about  22  Hertz  and  falls  off  to  values  one- 
third  the  peak  at  about  4  and  50  Hertz.  This 
spectrum  cube-root  scales  with  yield.  That 
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TABLE  1 


SUMMARY  OP  DURATION  ESTIMATES 


Bottom  SOFAR 


Model 

Depth 

D  (km) 

Range 

R(km) 

Duration 

(sec) 

Axis 

Velocity 

Extended 

Duration 

H 

5.0 

6600 

12 

1482.7 

55 

M 

5.0 

6600 

7 

1482.28 

32 

L 

5.0 

6600 

38 

1476.79 

55 

A 

o 

« 

in 

6600 

9 

1481.39 

38 

A 

4.5 

6600 

4 

1481.39 

33 

A 

5.5 

6600 

16 

1481.39 

45 

A  Mod 

in 

• 

6600 

7 

1481.39 

38 

A  Mod 

5.0 

6600 

12 

1481.39 

33 

L 

5.0 

3600 

22 

1476.79 

31 

A 

5.0 

3600 

6 

1481.39 

22 
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is,  at  2  KT  the  characteristic  frequencies  are 
a  factor  of  1.26  smaller.  The  spectral  esti¬ 
mates  do  not  account  for  the  recording  system 
response. 

•  The  duration  of  the  acoustic  signal  varies 
widely,  depending  on  the  assumptions  about 
the  travel  path.  Values  from  20  to  60 
seconds  seem  reasonable  for  a  range  of  about 
6600  kilometers. 

These  conclusions  are  based  on  an  analysis  that  includes 
many  approximations.  It  is  useful  to  list  these  here.  They 
are: 

•  The  explosion  is  represented  by  the  AFWL  1  KT 
nuclear  blast  standard.  The  near-surface  1  KT 
airblast  is  assumed  to  be  the  same  as  the 
airblast  from  2  KT  in  free  air. 

e  The  airblast  loading  is  applied  as  a  pres¬ 
sure  to  the  surface  of  the  ocean,  which  de¬ 
forms  elastically. 

e  The  direct  coupling  is  accounted  for  by  an 
equivalent  pressure  source  which  is  estimated 
via  the  process  described  in  Section  2.3. 

e  For  energy  travelling  along  paths  entirely 

within  the  ocean,  the  absorption  is  negligible 
at  the  low  frequencies  of  interest. 

•  The  travel  path  is  represented  by  a  single 
sound-velocity  profile.  Several  such  paths 
are  used. 

e  The  last  assumption  implicitly  assumes 
that  the  explosion  is  detonated  over  deep 
ocean. 

e  The  recording  system  response  is  not  in¬ 
cluded  in  our  spectral  estimates. 
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